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The intervortex interaction is investigated in very 
anisotropic layered superconductors in tilted magnetic field. 
In such a case, the crossing lattice of Abrikosov vortices (AVs) 
and Josephson vortices (JVs) appears. The interaction be- 
tween pancakes vortices (PVs), forming the AVs, and JVs 
produces the deformation of the AV line. It is demonstrated 
that in the result of this deformation a long range attraction 
between AVs is induced. This phenomenon is responsible for 
the dense vortex chains formation. The vortex structure in 
weak perpendicular magnetic field is the vortex chain phase, 
while in higher field a more complicated mixed vortex chain- 
vortex lattice phase emerges. 
PACS : 74.60.Ge. 

Vortex physics in layered superconductors occurred 
to be extremely rich and interesting. In moderately 
anisotropic superconductors, a tilted magnetic field leads 
to the formation of vortices inclined towards the super- 
conducting layers. The interaction between such tilted 
vortices happens to be quite unusual. 

In the plane defined by the vortex line direction and 
the c-axis (normal to the superconducting planes), the 
interaction between tilted vortices is attractive at long 
distances. Such attractive intervortex interaction is quite 
unexpected and leads to the formation of vortex chains, 
where the intervortex distance is governed by the bal- 
ance between the long range attraction and the short 
range repulsion. In layered superconductors, the exis- 
tence of these vortex chains in tilted field has been pre- 
dicted in [Ql and ||^, and subsequently confirmed by the 
decoration technique |^ , and the scanning-tunneling mi- 
croscopy [|j measurements, in YBa2Cu307 and NbSe2 
crystals respectively. The systems studied in [|[^] are 
characterized by a moderate anisotropy. The theoreti- 
cal approach ||,|^ is completely applicable to this case 
and gives a good qualitative and quantitative descrip- 
tion of this phenomenon. On the other hand, in the 
much more anisotropic Bi2Sr2CaCu202; (BSCCO) single 
crystals |^J^, a more complicated mixed vortex chain- 
vortex lattice phase has been observed. As it has been 
demonstrated 0, in strongly anisotropic layered super- 
conductors, in a tilted magnetic field, a crossing lattice 
of Abrikosov and Josephson vortices (JVs) , a "combined 
lattice" , must exist. The AV is in fact a line of pancake 
vortices (PVs) [p|-pd| interacting with JVs. Following 
Ip^ , the perpendicular vortex line formed by the PVs is 
deformed and attracted by JVs, so the JVs stacks ac- 
cumulate additional PVs, creating vortex row with en- 
hanced density |^J^ . This scenario has been proposed in 
p2[ to explain the mixed chain-lattice state formation. 



Very recently, using scanning Hall probe microscopy, the 
detailed studies of vortex chains in BSCCO have been 
performed and revealed the stability of the dense vortex 
chains state even in the absence of lattice and in very 
weak perpendicular magnetic field . 

In the present work, we demonstrate that the forma- 
tion of such dense vortex chain is due to the attraction 
between the deformed lines of pancakes vortices. The 
deformation, responsible for a long range attraction, ap- 
pears due to the interaction with JVs. In the result, the 
mechanism of vortex chains formation, in tilted field, in 
strongly anisotropic superconductors, appears to be quite 
similar to the case of moderately anisotropic supercon- 
ductors 

Further on, keeping in mind BSCCO, we will consider 
layered superconductors with high anisotropy ratio 7 — 
^c/^ab 200 — 500, where Ac is the penetration depth for 
currents along c-axis (perpendicular to the layers), and 
Xab is the penetration depth for currents in the ab plane 
(parallel to the layers) . The in-plane field — B cos 
penetrates inside the superconductor in the form of JVs, 
while the perpendicular field B^ = Bsind creates the 
PVs which interact with JVs via the Josephson coupling 
[|l2|,^. We consider the case of a very weak coupling 
of the layers when the Josephson's core radius Aj = 7s 
(s is the interlayer spacing), is lar ger than an in-plane 
penetration depth, i.e. Aj > Xab- In |l4l ; the limit of high 
field oriented near the (a, b) plane B ^ Hq = 00 / 
was considered and it was demonstrated that a zigzag 
displacement of PVs along x-axis is produced, see Fig. 
1. The authors of analyzed the case of the dense 
vortex lattice formed by PVs, but their method to treat 
the JVs and PV interaction permits us to calculate the 
shape of the single zigzag line of PVs too. 

In the limit of weak Josephson coupling, the interac- 
tion which stabilizes the straight PVs line is mainly of 
electromagnetic origin. Then, we may use the general 
expression for the energy of an arbitrary configuration of 
pancakes in the framework of the electromagnetic model 
B to calculate the energy increase due to the line defor- 
mation 
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where Ag/ / — -f^ , (k) is the Fourier transform of the 
total London vector of the n*^ layer $Ti(f) — ~ 
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) (the sum is over all the pancakes of the n' 
layer), $(k) = i(po , the function Gk — cosh (/cs) + 
^2\^^^fk ■ ^'-'^ zigzag deformation with an amplitude u, 
we have *2«(k) = e*"''*(k) and *2«+i(k) = e-'"i'*(k), 
and the energy increase per one layer occurs to be 
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On the other hand, the gain of the Josephson energy 



due to PVs displacements, calculated following 
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deformation energy 8Ef. 
its equilibrium value 



Minimization of the total 
5Ej with respect to u gives 
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and the limit ^ Hq corresponds to strongly overlap- 
ping Josephson cores. Note that the amplitude of the 
zigzag modulation satisfies the condition u <C X^^ , and 
then it is a relatively small deformation of the vortex line 
consisting of the PVs. 

Now let us calculate the interaction of such two zigzag 
vortex lines at the distance x, both located in the (x, z) 
plane, see Fig. 1. Those characteristic distances, which 
are of interest for us, are much smaller than Ac, and 
then in the calculation of the interaction energy we 
may neglect the Josephson coupling and use the gen- 
eral form of the energy in the pure electromagnetic 
limit (|l|). In our case, the London vectors of the lay- 
ers are *2n(k) = e™''*(k)(l + e"''''^) and *2n+i(k) = 
e^*"*'$(k)(l + e^*'''^). Performing in the necessary 
summation and using for Gk its expansion for ks <^ 1 
(which is perfectly justifiable when the distances of inter- 
est are larger than the interlayer distance s), we finally 
obtain the following expression for the interaction energy 
per one layer 
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Performing firstly the integration over ky and then over 
kx, we may present the interaction energy as 
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where Kq is the modified Bessel function of zero or- 
der. At long distances x ^ X^^, the Bessel function 

Kq decays exponentially, and the leading contri- 

bution comes from the last term in (^, which gives 



Eint{x) ~ — s (^ 2nX ) ' ^OTig distauccs the net 

interaction between the zigzag PVs lines is an attraction! 
At short distances, the interaction is repulsive, as usual 

Eint{x) ~ — s (^ 2^x ) (~^) ■ overall behavior of 
the interaction energy is presented in Fig. 2 for different 
deformation parameters e — ^ "^^"^ mini- 

mum energy may be easily found for small deformation 
parameter e using the asymptotic of the Bessel function 
Kq (z) « yj^exp{^z). With logarithmic accuracy the 
minimum realizes at 
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This result leads to the important physical conclusion 
that in the presence of Josephson vortices, the AVs, due 
to the long range attractive interaction, will form chains. 
The equilibrium distance between the vortices in the 
chain is given, with logarithmic accuracy, by the same 
expression (^ as in the two vortices case. In fact, the 
energy of the AV in chain is lower than the energy of 
a solitary vortex, and with the increase of the vor- 
tices will start to penetrate in form of chains. Further 
increase of B^ will increase the number of chains with- 
out changing in first approximation the distance between 
vortices in chains. Only when the distance between the 
chains reaches the value of the order of A^j, , the forma- 
tion of the usual Abrikosov lattice will occur. In the low 
field regime Bz ^ -ff^i = -r^ri- In ( -f^ ) , the well defined 



dense vortex chains must be present. 

Now let us discuss the influence of the parallel field on 
the vortex chain state. In the high field limit Bx 3> Hq, 
the increase of the parallel field decreases the deformation 
parameter, see (^, and thus, following (^, will slightly 
increase the intervortex distance inside the chain. A 
more important eflFect, actually, is the strong decrease 
of the potential dip in the intervortex interaction en- 
ergy (0) which will result in the melting of the vortex 
chain. The vortex chain occurs to be more stable at lower 
parallel field. However, the above approach is valid for 
Bx > Hq, and then the maximum stability of the chain 
is expected at Bx ~ Hq. At lower parallel field, the dis- 
tance between Josephson vortices along z-axis becomes 
much larger than 2s, as well as the distance between de- 
formed parts of the PVs line. Evidently, the attraction 
between AVs, due to its deformation, will weaken, and 
such a case needs a special analysis which is presented 
below. 

Let us consider the limit of weak parallel field, when 
the Josephson vortices are well separated and the dis- 
tance D between them along z-axis strongly exceeds the 
interlayer distance s, i.e. Bx ^ Hq, see Fig. 3. As it has 
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been demonstrated in |^2[, if the JV is located between 
the layers and 1, the pancake displacements u„ on the 



layer in the case of the single AV is 
2C„ 1 
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where Cn is a numerical coefficient of order one. Note 
that as it may be expected, this expression for n = 1 is 
of the same order of magnitude as (||) at the limit of its 
applicability at ~ Hq. 

To treat this situation, it is convenient to add a ficti- 
tious pair of pancake vortex and antivortex at the central 
line of the AV, see Fig. 3. Then, the obtained configura- 
tion will be equivalent to two straight Abrikosov vortices 
and to two vortex-antivortex pairs at the distance x in 
each layer. Vortex and antivortex are separated by a 
distance m„, and we are coming to the problem of the 
interaction of such magnetic dipoles. Firstly, note that 
the interaction between dipoles in the same layer is at- 
tractive, and it may be directly calculated with the help 
of (0): 
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It is much larger than the interaction between dipoles 
from different layers n and m, which may be attrac- 
tive or repulsive. From ([^) it may be estimated as 

\Eimi^)\ « 8^"^!^ fo"- s\n-m\ « A„,, i.e. 

ab afa 

containing the small additional factor ^ 1. More- 

over, for s |n — m| > A^^, this interaction decays expo- 
nentially. Taking into account that u„ varies as u„ ~ 
■7 — near the JV, we may finally conclude that the main 

contribution to the dipole interaction energy is coming 
from the interaction in the same layer, and it may esti- 
mated (per period D of the modulation of PVs line along 
z axis) as 
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The main contribution to the repulsion energy is coming 
from the straight AVs interaction, and per period D at 
distances x ^ A„^ is 
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As in the case of the dense JVs lattice, the prevailing in- 
teraction between AVs at long distances is an attraction. 



and it will lead to the vortex chain creation. However, 
as the attraction is coming only from the PVs near JVs, 
the relative strength of attraction is much smaller (other 
PVs contributing only to the repulsion). In the result, 
the distance between AVs in the chain may be estimated 
as 
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where e — x^V^ ~ iff ('^) ■ ^""^"^ energy gain 
due to the vortex chain formation per pancake is 
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i.e. it is decreasing with the decrease of B^- Then, in 
the case B^ <^ Hq, the perpendicular vortices also ap- 
pear as vortex chains, and when the perpendicular field 
increases it will lead firstly to the increase of the num- 
ber of these chains, each chain located at JVs. The dis- 
tance between the chains will be an integer number of 
the distance between JVs along y-axis. Finally, when 
all JVs will contain chains, with equilibrium distance 
Xmin between AVs, further increase of B^ will lead to 
the appearance of additional vortices in chains, decreas- 
ing the intervortex distance below Xmin- In such a case, 
the neighboring vortices repeal each other in the chain, 
but the repulsion energy is compensated by the gain of 
energy due to the trapping of the AVs by the JVs. How- 
ever with further increase of Bz , the repulsion energy will 
overcome the trapping energy and the formation of the 
usual Abrikosov lattice will start. Namely this case cor- 
responds to a mixed vortex chain- vortex lattice observed 
in 1^,^. It is interesting that in low perpendicular field 
Bz "C ^nd in the limit B^ < Hq the mixed vortex 
chain-vortex lattice may transform into a purely chain 
state with the increase of Bx- Indeed, the number of JVs 
increases, permitting therefore to accommodate all AVs. 
This type of behavior could be responsible for peculiar 
results on vortex lattice melting in the presence of JVs 
obtained in p^ . 

Comparing the results for low {Bx ^ Hq) and high 
{Bx 3> Hq) parallel field limits, we see that the distance 
between AVs in chain is always around A^,^ and slightly 
varies with B^- On the other hand, the energy of vortex 
coupling in a chain is maximal for Bx ^ Hq, and then 
they are more stable at this condition. At low or high 
Bx limits, we may expect the melting of vortex lines. 

Note also that the vortex chain state could also reveal 
an anisotropy of the critical current, the current along 
y axis will provoke the motion of the chain as a whole, 
while the current along x axis will provoke the vortex 
detachment from the chain. 

In conclusion, we predict a new qualitative effect, the 
long range attraction between AVs in the crossing vor- 
tex structure appearing in highly anisotropic layered su- 
perconductors. Such an effect is somewhat reminiscent 
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to the tilted vortex attraction in superconductors with 
moderate anisotropy If we consider the hne be- 

tween neighboring layers with tt phase difference as the 
center of a JV, the PVs displacement disrupts this line 
at the distance 2u. So we may consider the resulting 
structure as a tilted vortex line made from PVs and JVs 
parts (see the wavy line in Fig. 3). However the very 
important difference, with the straight tilted vortices, is 
that in the highly anisotropic case the vortex attraction 
is completely controlled by the parallel component of the 
magnetic field only. 
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of PVs stacks are at the distance x. JVs are presented 
by dashed lines. 

FIG. 2. The energy of interaction between two zigzag 
vortices as the function of the distance between them for 
different deformation parameter e = values. 

FIG. 3. Schematic picture of deformed AV lines in case 
of relatively weak parallel magnetic field {B^ ^ Hq) . 
The distance between the centers of AV lines is x. The 
JVs are presented by dashed lines. The period D of JVs 
lattice along the z-axis is much larger than the interlayer 
distance s. The PVs separation breaks the JV line, these 
parts are presented by cross, and the resulting "tilted 
vortex line" is depicted by a wavy line. On the upper 
layer, the fictitious vortex-antivortex pair is shown at the 
center of the AV. This procedure restores a straight AV 
line with additional vortex dipole presented as a shaded 
region. 
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Figure Captions 

FIG. 1. Zigzag deformation of the PVs stacks in high 
parallel magnetic field [B^ :$> Hq) , directed along the x- 
axis. u is the amplitude of deformation and the centers 



4 



A s{(l)jA7iX^Y 




